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ABSTRACT 

Small-scale fluctuations of magnetic field and velocity may be responsible for the Stokes asymmetries ob- 
served in all photospheric magnetic structures (the MISMA hypothesis). We support the hypothesis show- 
ing that atmospheres with optically-thin micro- structure reproduce the polarization of Fe I A6301.5 A and 
Fe I A6302.5 A observed in a sunspot. Ten thousand spectra were fitted by model MISMAs with two mag- 
netic components interleaved along the line-of- sight. Combining all the fits, we set up a semi-empirical model 
sunspot characterized by two components with very different magnetic field inclination. The major compo- 
nent, which contains most of the mass, is more vertical than the minor component. The field lines of the minor 
component are inclined below the horizontal plane throughout the penumbra. Magnetic field lines and mass 
flows are parallel, consequently, both upflows and downflows are present everywhere on the penumbra. Major 
and minor components have very different velocities (several hundred m s -1 for the major component versus 
10 km s -1 for the minor component), but the mass transported per unit time is similar. The similarity between 
the vertical mass flow and the magnetic flux of the two components suggests that field lines emerging as major 
component may return to the photosphere as minor component. If so, the observed magnetic field strength 
difference between components leads to a siphon flow whose magnitude and direction agree with the inferred 
Evershed flow. Several tests support the internal consistency of the retrieved model sunspot. The magnetic field 
vector B does not violate the VB = condition. The model sunspot reproduces the net circular polarization of 
the observed lines, plus the abnormal behavior of Fe I A 15648 A. The use of only one magnetic component to 
interpret the spectra leads to inferring upflows in the inner penumbra and downflows in the outer penumbra, in 
agreement with previous findings. 

Subject headings: line: profiles - polarization - methods: data analysis - Sun: magnetic fields - Sun: photo- 
sphere - sunspots 



1. INTRODUCTION 

Ubiquitous small-scale fluctuations of the magnetic field 
and the velocity may be responsible for the Stokes asymme- 
tries 1 observed in all photospheric magnetic structures. If 
the spatial scale is small enough to be optically-thin, one 
has a very efficient source of line as ymmetry. The existence 
of suc h gradients was conjectured by Sanchez Almeida et alJ 
(1996), that coined the acronym MISMA to describe 
magnetic atmospheres having optically-thin irregularities 
(MISMA = Micro-Structured Magnetic Atmosphere). If the 
conjecture turns out to be correct, we must understand the 
magnetic micro- structure for two fundamental reasons. First, 
many physical processes determining the magnetism at all 
scales take place at micro-scales (e.g., ohmic diffusion, ad- 
vection, onset of instabilities, viscous stresses, etc.). Second, 
it plays a key role for the proper interpretation of the mag- 
netic field measurements, which are prone to bias if the micro- 
structure is overlooke d. (For a more elaborated discussion on 
these two reasons, see Sanchez Almeida|2001[) 

This paper is the third of a series aimed at examining 
the conjecture in detail. We wanted to see (a) whether 
the observed asymmetries can be quantitatively explained by 
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1 We will use the standard notation for the Stokes parameters; / for the 
intensity, Q and U for the two independent types of linear polarization, and V 
for the circular polarization. The Stokes profiles are the function Stokes pa- 
rameter versus wavelength for a particular spectral line. These Stokes profiles 
are expected to be symmetric or anti- symmet r ic when the line is formed in a 
homogeneous atmosphere, e.g., Unno 1 1956); Landi Ded'Innocenti < 1992). 
Deviations from these symmetries are called Stokes asymmetries. 



MISMAs, and (b) what are the physical properties of the 
micro- structure. The first paper introduces an inversion code 
(IC) based on the MISMA concept. It is a (fairly) auto- 
matic algorithm to fit Stokes profiles using synthetic profiles 
formed in atmospheres with optically-thin micro- structure 
( Sanchez Almeida 1997, Paper I). As a test case, Pa- 
per I reproduces the asymmetries found in plage regions by 
Stenfloetal. (1984). The second paper describes the suc- 
cessful application of the IC to quiet Sun network and inter- 
network r egions (|Sanchez A lmeida & Lites 2000, Paper II ; 
see also Socas-Navarro & Sanchez Almeida 2002). The 
present paper extends the MISMA inversions to a sunspot, 
showing that the micro- structure is also able to cope with the 
Stokes asymmetries of sunspots. 

The paper is organized as follows. The properties of the 
spectro-polarimetric data subject to analysis are summarized 
in § 13 Section |3] describes the constraints used to set up 
the model MISMAs, the quality of the fits, and the proce- 
dure to assemble independent inversions from different pix- 
els into a single coherent structure. The main results are set 
forth in the seven subsections of § |4| Specific aspects of the 
resulting model sunspot are examined in § |5] (magnetic flux 
conservation) and §|S|(mass conservation). Section § CI ex- 
plains the physical origin of the asymmetries from a radiative 
transfer point of view. Finally, § comments on and sum- 
marizes the contents of the paper. The interpretation of these 
results in the context of the global penum bral magnetic field 
topolo gy is deferred for a separate paper ( Sanchez Almeida 
2005). We have preferred to disconnect the fact that MISMAs 
account for the Stokes asymmetries in sunspots (this work), 
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Fig. 1 . — Continuum image of the inverted sunspot. Coordinates are given 
in Mm from the adopted sunspot center. Those pixels which have been in- 
verted are artificially enhanced so that the sunspot and its surroundings stand 
out against a dark background of pixels without inversion. The white contour 
shows the border of the sunspot, and it is used for reference along the paper. 



from the more speculative physical interpretation of the re- 
trieved micro- structure. Even if the physical interpretation is 
eventually proved to be wrong, the micro- structure would still 
account for the asymmetries. 

2. OBSERVATIONS AND DATA 

The Stokes /, Q, U and V profi les analyzed in this work 
were obtained bv lLites et alJ (119981) wi th the Advanced Stokes 
Polarimeter (ASP; Elmore et al. 1992). They correspond to 
the sunspot in active region NOAA 7912 (Fig.[IJ observed on 
1995 October 15 when passing close to the solar disk cen- 
ter (S5, E9, cosine of heliocentric angle \i = 0.96). The re- 
gion was scanned with a step-size of 0."75 in the East- West 
direction, perpendicular to the slit. The pixel size along the 
spectrograph slit 0." 37, is enough to sample a seeing close 
to (but better than) 1". The Field-Of-View (FOV) is set by 
the length of the slit, and the size of the scan, and it turns 
out to be 85" x 135". The ASP provides the polarization of 
Fe I A6301.5 A and Fe I A6302.5 A every 12.8 mA. An inte- 
gration time per slit position of 10 sec affords a noise per spec- 
tral pixel better than 10~ 3 times the continuum intensity. The 
spectra are corre cted for instrumental po larization to a level 
below this noise JSkumanic h et aljfl997l) . For furt her details 
on the observational setup and the target, refer to iLites et all 
( 1998). We use this rather old data set because its image qual- 
ity was among the best ASP maps available when the work 
began in 1999. 

The sunspot has negative polarity, i.e., the magnetic field 
lines are anti-parallel to the vertical direction at the umbral 
core. This fact introduces some minor complications in the 
interpretation of the sunspot magnetic field. They will be 
pointed out and clarified whenever they interfere with our ar- 
gumentation. 

The zero of the absolute velocity scale has been set using 
the Stokes V zero-crossing wavelengths in the plage around 
the sunspot. We use its mean position as null velocity. The 
same procedure was used in Paper II since observations have 
shown that the zero-crossing wavelengths of Stokes V pro- 
files in pl age are very close to the laboratory wave lengths of 
the lines (Solanki 1986: lMartinez Pillet et aj1T9 97). The sys- 
tematic error of our absolute velocities should be less than 
0.25 km s" 1 . 
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Fig. 2. — Top: cartoon with the kind of atmosphere used for inversion of 
sunspot Stokes profiles. It represents magnetic field lines in one resolution el- 
ement. They overlap along any line-of-sight . (a) A view from above showing 
how the field lines are preferentially aligned along the sunspot radial direc- 
tion, (b) A view from the side showing field lines with very different vertical 
inclinations. 

3. INVERSION 

3.1. MISMA scenario to represent umbrae and penumbrae 

We employ the scenario proposed in 

Sanchez Almeida et al ] (11993 § 3.1), which resem- 
bles a scaled-down versio n of the uncombed fields of 
Solanki & Montavon (1993). The atmosphere has two 
magnetic components which differ in velocity, inclination, as 
well as in many other physical parameters. It is a MISMA, 
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therefore, each line-of-sight crosses several different sub- 
structures (see the sche matic in Fig.|3. Since we will use the 
Inversion Code (IC) by ISan chez Almeidaj (119971) . these two 
components meet several physical requirements. 

1. The temperature of the two magnetic components is 
the same at every geometrical height. We adopt 
this approach for consistency, since the radiative ex- 
change among the optically-thin micro- structures of the 
MISMA s hould efficiently smooth t emperatures in sec- 
onds (e.g., Sni egelH 957tlstixll1 991h 2 . 

2. In addition to the two magnetic components, there is 
a third non-magnetic component. It is introduced for 
techni cal r easons and plays no central role. As we will 
see (§ 14.1b . the fraction of photospheric volume that it 
occupies within the sunspot is always marginal. 

3. The three components are in lateral pressure balance at 
every geometrical height (i.e., they have the same to- 
tal pressure, which is the sum of the gas pressure plus 
the magnetic pressure). This constraint and the assump- 
tion of hydrostatic equilibrium along field lines provide 
the full vertical stratification of magnetic field strengths 
from the strengths at one height. The field strengths of 
the different components are different. 

4. The two magnetic components have different inclina- 
tions with respect to the vertical direction, but the same 
azimuth. Inclinations and azimuths do not vary with 
height in the atmosphere. 

5. Only motions along magnetic field lines are allowed. 
The mass conservation along field lines determines the 
velocity stratification by setting the velocity at a given 
height. 

6. Unpolarized stray light is accounted for. Dealing with 
stray light introduces three free parameters: the frac- 
tion of stray light, a global velocity, and a microturbu- 
lence. The thermodynamic stratification is that of the 
non-magnetic component of the MISMA model atmo- 
sphere. 

The actual model atmospheres are described using 23 pa- 
rameters: the temperature of the non-magnetic component at 
four heights, a global scaling factor for the temperature of 
the magnetic components (item [I] above), two magnetic field 
strengths at the bottom of the atmosphere (or BID, which 
stands for Bottom of Integration Domain) (item|3]>, two oc- 
cupation fractions at BID (i.e., the fraction of atmospheric 
volume occupied by each magnetic component; item 01, two 
velocities at BID (item |3, two magnetic field inclinations 
(item|4j>, one magnetic field azimuth (item|4|>, four parame- 
ters to define the distribution of microturbulent velocities (two 
for the magnetic components, and two for the non-magnetic 
component), three parameters to account for the stray light 
(item|6|>, one macro-turbulent velocity (which simulates in- 
strumental broadening plus macro-turbulence of solar origin) 
and, finally, the maximum pressure at BID. The latter varies 
to assure that the model atmosphere embraces the region be- 
tween logr c = 1 and logr c = -4 (with r c the continuum opti- 
cal depth at 6302 A). This range of optical depths guarantees 

2 However, this hypothesis makes our model insensitive to the temperature 
variations associated with the smallest penumbral filaments observed by, e.g., 
Schar mereTaTI i200l) . 



that the lines Fe I A6302.5 A and Fe I A6301.5 A are entirely 
formed within the bounds of the model atmosphere. Recov- 
ery of the complete at mospheric structure from t hese 23 pa- 
rameters is described in Sanchez Almeida ( 1997) 3 . These 23 
parameters must be set during the least- squares fitting proce- 
dure. Although the number of free parameters seems large, it 
has to be compared with the 536 different observed quantities 
to be fitted (67 wavelengths per line and per Stokes parameter; 
see Fig.0. These observables are not fully independent (e.g., 
the finite spectral resolution links neighbor wavelengths, and 
the radiative transfer couples the four Stokes parameters), but 
the author believes that the number of observables is so large 
compared with the number of free parameters that the number 
of independent observables almost certainly exceeds 23. 

Figure|2|illustrates the kind of complex distribution of mag- 
netic field lines that our schematic model MISMAs may be 
modeling. Field lines like these are represented by our IC as, 
e.g., the model MISMA in Fig.|4| If the micro- structure of 
the atmosphere is optically-thin then the emerging spectrum 
does not depend o n details of the atmosphere bu t only on its 
average properties ( Sanchez Almeida et al. 1996). Therefore, 
describing a messy bunch of field lines with only two compo- 
nents is consistent with our MISMA assumption. Obviously, 
we cannot infer details of the micro- structure but volume av- 
erage properties. In particular, we are insensitive to the size 
of the micro- structures, which we consider to be optically- 
thin and so smaller than 100 km. The cartoon in Fig. |2| rep- 
resents magnetic field lines existing in a single resolution ele- 
ment. The field lines are aligned approximately radially from 
the sunspot center, having different azimuths and inclinations. 
Figure^ shows a view from the top of the 3 -dimensional dis- 
tribution of field lines. They show a global radial orientation 
but, due to slight misalignments and swings, many different 
field lines overlap along any line-of-sight. The axes, given in 
Mm, are for guidance only. Figure^ represents a side view 
of the same distribution of field lines. Different magnetic field 
inclinations coexist in the resolution elements. 

The MISMA scenario and the requirements 
used for sunspot inversions are basically those 
that successfully reproduce the Stokes profiles ob- 
served in plage , network and internetwork r egions 
JSanchez Almeida H 9971: ISanchez A lmeida & Litesl [2000; 
Socas-Navarro & Sanchez Almeida 2002). The main differ- 
ence being the inclination of the magnetic fields, and the 
different temperatures of the magnetic components and the 
background. The similarities between scenarios probably 
reflects both the degeneracy of the radiative transfer for 
MISMAs and the fact that the employed MISMA scenario 
grasps the essential degrees of freedom required to reproduce 
most asymmetries found in the photosphere. 

As explained above, there are two magnetic component in 
each model MISMA. They will be denoted as minor and ma- 
jor depending on which one has largest ap, i.e., the largest 
product of occupation fraction a times density p. Conse- 
quently, the major component is the one containing most of 
the mass of the atmosphere. This intuitive classification is 
also used in Paper II. Here, however, the component of largest 
ap may change along the atmosphere, making the definition 
ambiguous. Fortunately, the change occurs in a small fraction 
of inversions (3%) and, even in these cases, it only affects a 
few superficial layers. (The classification is based on ap at 
the bottom of the atmosphere.) 

3 But note the misprint in equation (13). 
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Fig. 3 . — Examples of the quality of the fits provided by the MISMA inversions. They show Stokes /, Q, U and V profiles of Fe I A6301 .5 A and Fe I A6302.5 A. 
(The different types of profiles are separated by vertical solid lines, and the division between Fe I A6301.5 A and Fe I A6302.5 A is indicated by vertical dashed 
lines.) The solid lines are synthetic profiles whereas the dotted lines show observations, (a) Profiles in the limb-side penumbra, (b) Profiles in the umbra, (c) 
Profiles in the center-side penumbra, (d) Profiles outside the limb-side penumbra. The abscissae represent the pixel number in an array used to carry out the fits. 
Within each type of Stokes profile, wavelengths increase to the right, with a gap between Fe I A6301.5 A and Fe I A6302.5 A to exclude a telluric line. Note the 
large asymmetries of some of the profiles, and the fair fit provided by the inversions; e.g., Stokes V in (a) and (d). 



3.2. Strategy for inversion 

First, the center of the sunspot was set by trial and error. 
We minimize the departures from axi- symmetry of the distri- 
bution of observed continuum intensities. The chosen center 
is 1.3 Mm away from the darkest point in the umbra. The re- 
sulting intensity distribution is fairly symmetric, leading to a 
sunspot radius of about 17 Mm, which is the value adopted in 
the paper. 

We choose for inversion all those pixels within 24 Mm from 
the sunspot center having a degree of polarization larger than 
0.5% (average polarization over a bandpass 75 mA wide in 
the blue wing of Fe I A6302.5 A). This criterion yields 10130 
different pixels encompassing the entire sunspot (Fig.[IJ. We 
have carried out an inversion for each one of these pixels u sing 
the IC described in Paper I and the scenario presented in $ 13.11 
The strategy for inversion is slightly different from the one in 
Paper II. We also begin by classifying all Stokes profiles us- 
ing a principal component analysis scheme. This procedure 
sorts the observed Stokes profiles according to similarities in 
shape. We start the inversion at the sunspot center, and then 
proceed with the closest profile according the classification. 
However, rather than using the last retrieved atmosphere to 



initialize each new inversion, we look up in a data base with 
all previously inverted synthetic profiles. The atmosphere of 
the closest one (in a least squares sense) is chosen for initial- 
ization. The process is repeated to completion of the 10130 
sets of Stokes profiles. Each automatic inversion was visually 
inspected before acceptance. Failures were re- started from a 
different model atmosphere until the agreement was satisfac- 
tory. 

3.3. Quality of the fits 

Figure [3] shows four examples of the quality of the fits 
provided by the MISMA inversions. We have chosen the 
Stokes profiles at four positions on the sunspot, (a) limb-side 
penumbra, (b) umbra, (c) center- side penumbra, and (d) out- 
side the penumbral border defined by the continuum image 
(the letters correspond to the labels on the figure, whereas 
limb- side and center- side refer to the portion of penumbra 
closest to or farthest from the solar limb). These four re- 
gions produce Stokes profiles with the variety of shapes to 
be reproduced with the inversions, including severe Stokes 
asymmetries. The Stokes V profile in Fig Eh shows the 
so-ca lled cross-over effect (Grigorie v & Katzffl972tlGo lovko 
with Stokes V having three or more lobes. It occurs 
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Fig. 4. — Typical model MISMA retrieved by the IC used in this work. It 
reproduces the profiles in Fig. [3^. The different panels represent the varia- 
tion with height in the atmosphere of various physical parameters (from left 
to right and from top to bottom): the continuum optical depth, the magnetic 
field strength, the temperature, the gas pressure, the density, the electron pres- 
sure, the occupation fraction, and the speed along field lines. Except for the 
upper-left panel, the type of line identifies the three components of the model 
MISMA; the solid lines correspond to the non-magnetic component and the 
dotted and dashed lines to the magnetic components. The diamonds indicate 
the free parameters of the model, from which the entire stratification of the 
atmosphere has been deduced. Note that the atmosphere is only meaningful 
within the range of optical depths where the observed lines are formed (e.g., 
the drop of major component magnetic field strength below log r c = is not 
significant; logr c = occurs at 100 km). 



in the neutral line, where the mean magnetic field is per- 
pendicular to the line-of- sight, and Stokes V should be zero 
if the magnetic field were spat ially resolved (Makita 1986; 
Sanchez Almeida & Litesl fT992). Despite the quality of the 
fits, we would like to point out that they are not as good as for 
the quiet Sun network and internetwork profiles describ ed in 
Paper II and lSocas-Navarro & Sanchez Almeidal (|2002). We 
have not been able to pin down the cause, however, it ob- 
viously implies that the scenario adopted for inversion suits 
better the conditions in the quiet Sun. 

Yet another hint of the quality of the fits is provided by the 
ability to reproduce the observed net circular polarization, a 



sunspot axis 
(vertical direction) 




Fig. 5. — Schematic showing the two possible solutions compatible with a 
given Zeeman absorption. The two have the same magnetic field component 
along the line-of-sight, and differ by 180° in azimuth. We have tried to select 
the solution that converges toward the sunspot axis (solution 1). (We deal 
with a negative polarity sunspot, so that solution 1 tends to be anti-parallel to 
the vertical direction.) 



result analyzed in detail in § 14.61 

3.4. Removal of the azimuthal ambiguity 

A rotation of the magnetic field vector by 1 80° about the 
line-of-sigh t does not modify the Zeem an absorption (e.g., 
Unno 1956; Landi Deel'Innocenti 1992). This leads to the 
well known azimuthal ambiguity of the magnetic field de- 
rived from Zeeman magnetometry, which has to sorted out 
before proceeding. Since the MI SMA model atmospheres 
have two magnetic components (§ 13.1b . each one with this 
1 80° ambiguity, each pixel has four different magnetic con- 
figurations compatible with the observation 4 . Among them, 
we chose the one that best fits in a global magnetic field that 
converges toward the sunspot core, and is as smooth as possi- 
ble. First, we begin with the major components of all pixels 
( i.e., with the components containing most of the mass; see 
33.11) . The smoothest solution is obtained after an iterative 
procedure where the direction of the magnetic field vector in 
a given pixel is compared with the direction of the mean field 
provided by the surrounding pixels. Then the solution closest 
to the mean direction is chosen. This procedure is applied to 
each pixel and then iterated to a point where no pixel over the 
full FOV has to be modified. The iterative procedure starts by 
choosing that solution representing a magnetic field converg- 
ing toward the sunspot core (see Fig.|3>. If the two solutions 
of a pixel are of this kind, then the closest to the sunspot radial 
direction is chosen. We tried various box sizes to compute the 
local average (from 3x3 pixels to 15x15 pixels). The disam- 
biguation of the sunspot magnetic field is almost insensitive to 
this parameter, which mostly affects the field outside the visi- 
ble sunspot. We finally adopted 10x10 pixels, which provides 
a fair compromise. Once the map of the major component is 
retrieved, we choose that minor component solution whose 
direction is closest to the major component direction. The al- 
gorithm sketched above provides a fairly axi- symmetric mag- 

4 Two remarks are in order. First, neither the azimuth nor the inclination 
of each component vary with height in the atmosphere ($ 13.11 . Having vary- 
ing azimuths and inclinations would make the ambiguity intractable. Sec- 
ond, the Doppler shifts only depend on the projection of the velocity along 
the line-of-sight. Since we assume magnetic-field-aligned motions, all these 
four solutions have the same velocity along the line-of-sight and therefore are 
compatible with the inferred velocity field. 
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netic field distribution centered at the sunspot core inferred 
from the continuum intensity. Figure |6| shows maps of incli- 
nation and azimuth for the major components. Figured shows 
the same maps for the minor components. 

3.5. Setting a common scale of heights for the full inversion 

Each pixel has been inverted independently from the rest. 
Therefore each one has its own scale of heights. Part of the 
analysis of the forthcoming sections becomes more meaning- 
ful when carried out in a common scale of heights for the 
whole sunspot. We have set such common scale assuming 
that the total pressure (gas pressure plus magnetic pressure) 
is similar among neighbor pixels. (This approach was suc- 
cessfully applied in Paper II ; § 4.2.) The absolute vertical 
scale of each pixel is set by comparison with the pressure of 
the nearest neighbors whose absolute scale is already known. 
We minimize the differences between the logarithms of the 
total pressures allowing only for a global shift of the scale of 
heights. The procedure starts at the center of the sunspot and 
then continues to the penumbral border, and beyond. The zero 
of the absolute scale is set by the layer r c = 1 in the stray light 
component of the model atmospheres far from the sunspot 
(further than 23 Mm; see Fig.[l}. Theses model atmospheres 
represent the non-magnetic Sun so that the absolute heights 
are referred to the base of the non-magnetic photosphere. The 
uncertainty introduced by the assumption that the total pres- 
sure does not change between adjacent pixels is discussed in 
Appendix lAl It is smaller than 100 km for determining the 
Wilson depression of the umbra, and it becomes negligible 
small when computing pixel-to-pixel variations (< 4 km). 

4. MAIN RESULTS OF THE MISMA INVERSION OF A 
SUNSPOT 

4.1. Occupation fraction 

As described in § 13.11 the model MISMAs used for inver- 
sion have three different components, two magnetic plus a 
non-magnetic background 5 . The occupation fraction (OCF) 
is the fraction of atmospheric volume occupied by each one 
of them. In principle the OCFs vary with height in the at- 
mosphere (Paper I) . However, height variations are expected 
only when the OCF of the background differs from zero, since 
the magnetic OCFs grow at the expense of non-magnetic vol- 
ume. Since the OCF of the background is very small (see 
below), the OCF of the magnetic components remains basi- 
cally constant with height. Figure[S] shows the retrieved OCFs 
as a function of the radial distance from the sunspot center. 
We plot the values at the height in the atmosphere where 
logr c = -1, an optical depth where most of the atmospheric 
information extracted from the inversion comes from. The 
shaded areas represent plus and minus one standard deviation 
among all the points at the same radial distance. Note that the 
major component has the largest OCF all the way from the 
umbra to the sunspot surroundings. This result is not a triv- 
ial consequence of the definition of major a nd m inor, since 
the definition also involves mass densities (§ 13.1b . The ma- 
jor component OCF goes from more than 90% in the umbra 
to about 70% in the outer penumbra (the solid line in Fig.©. 

5 The non-magnetic backgroun d sho uld not be confused with the stray light 
component included in the fits (§ 13.11 item ^|6|. The stray light component 
will not be discussed in the paper since it shows a rather expected behavior 
similar to the non-magnetic stray light component deduced from the Milne- 
Eddington inversions of the same sunspot. The stray light fill factor goes 
from 2 % in the umbral core to 25 % at the penumbral border. 



The minor component is significantly different from zero only 
outside the umbra (from 5 Mm on; Fig. [8] the dashed line). It 
occupies 30% of the penumbral volume. 

4.2. Magnetic field inclinations and azimuths 

Figure[6]shows inclinations and azimuths of the major com- 
ponent. Figure shows maps of the same parameters for the 
minor component. The actual values are better appreciated 
in Fig. |9| where we represent the mean azimuths and incli- 
nations, as well as the standard deviation of these quantities 
among those points at the same distance from the sunspot cen- 
ter. (For a visual definition of vertical inclination and relative 
azimuth, see Fig. [K)|) Both the major component and (to less 
extent) the minor component follow the radial direction. The 
relative azimuth with respect to the radial direction is of the 
order of 180°, as it should be for a radial magnetic field in a 
negative polarity sunspot (Fig. |9j)). (The properties of the mi- 
nor component for distances smaller than 5 Mm are excluded. 
The minor component of these points does not contribute to 
the observed signals and therefore, any inference based on it 
results unreliable.) The magnetic field inclination of the ma- 
jor component follows a variation with radial position typical 
of sunspots (from say, 170° at the sunspot core, to about 15° 
of f the horizontal a t the sunspot boundary; see., e.g., Fig. 9 
m iLites etal.lfT99l . Together with this smooth and expected 
behavior, the minor component is (with large scatter but be- 
yond any doubt) always less inclined than the horizontal (see 
the dashed line in Fig. |9t). We deal with a negative polar- 
ity sunspot and, therefore, visualizing the magnetic topology 
from the inclination is not so straightforward. The cartoon 
in Fig. [H] has been laid out to illustrate the magnetic topol- 
ogy. This inclination of the minor component is common to 
the whole penumbra, implying that magnetic field lines are re- 
turning to the solar interior over the entire penumbra. This is 
possibly the biggest surprise of these MISMA invers ions . As 
we will discuss in connection with the mass flows ($ 14. 51 . in- 
ferences of magnetic fields returning to the sub-photosphere 
at the outer penumbra and beyond do exist i n the literature 
(Rimmele 1995; Westendorp Plaza et al. 1997). The presence 
of magnetic field lines that return throughout the penumbra is 
new, though. 

4.3. Magnetic field strengths and fluxes 

Figure[l2| shows the variation of the magnetic field strength 
with the distance from the sunspot center. The field strength 
varies with height, so we plot the value at logr c = -1. The 
two magnetic components are shown separately. (The field 
strength of the minor component in the umbra is not included 
since its occup ation fraction is negligible and therefore its 
value dubious; $ 14.11 ) The field strengths go from 3 kG at the 
sunspot core to 1 kG at the transition between penumbra and 
its surroundin gs. The observed v ariation is typical of sunspots 
(e.g., Fig. 9 in lLites et al.lfl9 93). More interesting is the sys- 
tematic difference of field strength between components. The 
mean field strength of the minor component is about 200 G 
larger than that of the major component. This excess is signifi- 
cant. The dispersion among the field strengths of all the pixels 
at a radial distance is also of the order of 200 G (the shaded 
regions in Fig.[T2l>. The dispersion of the mean is reduced by 
a factor square root of the number of pixels contributing to 
each radial point. Typically, several hundred pixels contribute 
to the mean so that the uncertainty of the mean is only a few 
tens of G, and so, well below the observed 200 G difference. 
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Fig. 6. — Maps of the inclination (left) and the local azimuth (right) of the major magnetic component of the MISMA inversions. The inclinations are scaled 
between 0° and 180° whereas the azimuths go from 0° to 360°. Note the axi-symmetry of the retrieved magnetic fields. The white outline corresponds to the 
external penumbral border (see Fig.[l}. 
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Fig. 8. — Radial variation of the fraction of resolution element occupied 
by each component of the model MISMAs. The radial distance is given in 
Mm, and all occupation fractions are evaluated at logr c = -1. The shaded 
regions represent the standard deviation among the points at a distance from 
the sunspot core. 



The error bars of the mean are also included in Fig. [12] but 
they go almost unnoticed due to their small values. 

The excess of field strength creates a surplus of magnetic 
pressure and, consequently, a deficit of gas pressure in the mi- 



nor component. Let us assume that the minor and major com- 
ponents of the same or neighbor resolution elements are con- 
nected by field lines. (We have not been able to come up with 
any other possibility to understand the presence of field lines 
pointing up and down throughout the penumbra. However, 
this connection between the field lines of the minor and major 
components does not result from the inversions and is specu- 
lative.) Then the gradient of gas pressure between minor and 
major components is bound to drive a mass flow. This is noth- 
ing but the siphon flow mechanism invoked to explain the Ev- 
ershed flow (iMever & Schmidtlll968t iThomas & Montesinosl 
The seemingly small difference of field strength that 
we infer can drive a very fast flow, as attested by a simple 
order of magnitude estimate. For a typical magnetic field 
strength in the penumbra B of 2 kG, 



AB/B^0.1. 



(1) 



The square velocity difference produced by a gas pressure ex- 



cess AP g is 



AU 2 



2APn 



(2) 



where we have use d the equatio n of motion for a horizontal 
thin flux tube (e.g.. ISpruitll 19811) . Since the gas pressure dif- 
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Fig. 9. — Radial variation of the magnetic field inclination (a) and relative 
azimuth (b) of the two magnetic components of the model MISMAs. The 
lines represent the mean values among those pixels at a distance from the 
sunspot center. The shaded regions correspond to the standard deviation. 
Each type of line identifies one component, as described by the inset in (a). 
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Fig. 11. — Schematic showing the direction of the major and minor com- 
ponents relative to the vertical and radial directions. It shows the components 
at three different radial distances. Since we deal with a negative polarity 
sunspot, the magnetic field points downward at the sunspot axis. 
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Fig. 12. — Radial variation of the magnetic field strength of the two mag- 
netic components. The mean field strength of the minor component is 10% 
larger than the strength of the major component. The shaded regions repre- 
sent the standard deviation among all the points at a distance from the sunspot 
core. The lines correspond to the mean values. Error bars for these mean val- 
ues are also included, but they go unnoticed due to their small values. 
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Fig. 10. — Schematic to show the inclination and relative azimuth used to 
present the results of the inversion. It displays the coordinate system used to 
describe the magnetic field vector (the thick arrow) at a given radial distance 
from the sunspot core. The relative azimuth is the azimuth in a local cylindri- 
cal coordinate system where the radial direction is the radial direction from 
the sunspot center. The inclination is inclination with respect to the vertical 
direction. The coordinate system is different for different pixels. For the sake 
of clarity, we show a magnetic field of positive polarity, however, the true 
sunspot is of negative polarity. For this reason the relative azimuths tend to 
be 180° (Fig.0>). 



ference is due to the magnetic field strength difference, then 

B 2 AB 

AP„ = . 

8 4tt B 

Using equations and one finds, 

B 2 Afl-iV2 



(AU l ) 



2.1/2 ( 



2itp B 



(3) 



(4) 



with the symbol p standing for the dens ity. F or typical densi- 
ties at logr c ~ -1 (~ 10~ 7 g cm -3 ; see $ 14.41) . 

(A*/ 2 ) 1/2 ~8kms _1 . (5) 

This velocity is similar to the m inor component velocities in- 
ferred from the inversions; see $ 14.51 

According to the current paradigm, a large fraction of 
the sunspot magne tic flux emerges in the penumbra (e.g. 
Solanki & Schmidt 1993). It leads to the concept of deep 
penumbra, where the boundary between magnetized plas- 
mas and non-magnetic environment occurs well below the 
observed photospheric layers. To our surprise, the model 
sunspot turns out to have very little vertical magnetic flux 
across the penumbra. Our finding does not imply a shallow 
penumbra, though. It follows from the large amount of return 
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flux in the penumbra. The two magnetic components tend 
to hav e vertical magnetic field strengths with opposite signs 
(§ 14. 21 and Figs. 13 and [HJ. This leads to a cancellation of the 
average vertical component of the field strength <B Z >. (The 
angle brackets < > are used in the paper to represent vol- 
ume average, i.e., an average weighted with the OCFs.) The 
mean vertical component at logr c = -1 is shown in Fig.[T3k. 
It tends to zero in the external penumbra (say, for distances 
between 10 Mm and 16 Mm). However, the unsigned mean 
< \B Z \ > is certainly different from zero. Such difference be- 
tween signed mean and unsigned mean produces significant 
differences between the signed and unsigned magnetic fluxes 
of the sunspot. The signed vertical magnetic flux of those 
points with radii between r\ and r 2 is 



$ z (ri <r<r 2 ) = 2ir r< B z >dr, 

J n 
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(6) 
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with <B Z > the azimuthal average of < B z >, 



<B Z > dx, 




(7) 



and x the absolute azimuth of each point. (The unsigned mean 
is formally identical to equation (§] replacing < B z > with 
< \B Z \ >.) FigurefT3b shows the integrand of equation ®, i.e., 
the contribution of each radius to the sunspot flux. It is clear 
how the penumbra makes a significant contribution to the un- 
signed flux of the sunspot, but little to the signed flux. Table[T] 
quantifies how much. It lists the vertical magnetic flux of the 
sunspot, the umbra and the penumbra. The signed flux of the 
sunspot is 33% smaller than the unsigned flux. Flux cancel- 
lation primarily occurs in the penumbra (58% reduction) and 
very little in the umbra (9% reduction). The cancellation at the 
penumbral border is more than 80%, a result which could be 
naturally explained if field lines emerge, bend over, and return 
in a distance similar to the angular resolution. (A detailled 
discussion about the physical interpretatio n of such specula- 
tive a rguments is deferred to another paper; Sanchez Almeida 
l2005h The penumbra has as much unsigned flux as the um- 
bra, but only 40% of the signed flux. Table[T]also includes the 
fluxes of a standard model sunspot without return flux. The 
flux of the standard model sunspot lies in between our signed 
and unsigned fluxes. 

4.4. Mass Density 

Figure[l4|shows the variation of the density along the radial 
direction at logr c = -1. The densities of the major component 
are not very different from the densities of the quiet Sun at 
this optical depth. Note, however, how the density of the mi- 
nor component is five times smaller. The association of minor 
component with low density is partly due to the definition of 
minor and major, which involves the density. In addition, the 
enhanced magnetic field strength of the minor component de- 
scribed above demands a decrease of gas pressure and, con- 
sequently, a decrease of density. (Recall t hat m ajor and minor 
temperatures are forced to be identical; $ 13.11 ) 

4.5. Motions and the Evershed effect 

All motions are aligned along the magnetic field lines. 
T his is a constraint directly imposed on the model MISMAs 

(§ED- 

The velocity pattern retrieved from the inversion is simple 
and regular; see Fig. [21 The minor component, with the 
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Fig. 13. — (a) Volume average of the vertical magnetic field strength. The 
solid line shows the signed mean whereas the dashed line corresponds to the 
unsigned mean. Note the large difference in the external penumbra (say, for 
distances between 10 and 16 Mm.) The shaded regions represent plus-minus 
one standard deviation among all those points at the same distance from the 
sunspot center, (b) Average vertical field times distance to the sunspot center 
(r = Distance). The integral of this quantity between two radii provides the 
vertical magnetic flux. 
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Fig. 14. — Radial variation of the density of the two magnetic compo- 
nents. The major component is five times denser than the minor component. 
Distances are given in Mm from the sunspot center. 



lowest mass density and the highest magnetic field strength, 
has very large velocities. They typically reach 10 km s -1 at 
logr c = -1. The major component shows much smaller val- 
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TABLE 1 

Vertical magnetic flux across the sunspot 









Unsigned 


Signed 


Reference a 


<S> Z (0 <r<R s ) 


full sunspot 


$ 


= 7.3 x 10 21 Mx 


0.67 $o 


0.99 $o 


$z(R S /2 <r<R s ) 


penumbra b 


$1 


= 0.49 <£> 


0.42 $i 


0.83 $i 


<£ Z (R S /V2 <r<R s ) 


outside penumbra c 


<£> 2 


= 0.25 $o 


0.20 $ 2 


0.63 $ 2 


$ z (0 < r < R s /2) 


umbra 


$ 3 


= 0.51 3> 


0.91 3> 3 


1.14 $ 3 



a Normalized sunspot model in Solanki (2000), Table 14.34, with a maximum vertical field 
strength of 2.8 kG and the radius of our sunspot. 
b We assume the umbral radius to be half of the sunspot radius Rs- 
c Penumbral rim having half of the sunspot area. 



ues. Both velocities are directed radially outward along the 
field lines. Since the magnetic field lines are not horizontal 
(see, Fig.[9|and § I4.2> the velocities have both horizontal and 
vertical components. Figure [16| shows the radial velocity U r 
and vertical velocity U z of the major (top) and minor (bottom) 
magnetic components at logr c = -1. The vertical velocity of 
the major component is positive, meaning that it is directed 
upward. The vertical velocity of the minor component is di- 
rected downward. This downward directed velocity occurs all 
over the penumbra (Fig. [EJ. The small upward motions in 
the umbra (Fig. [l6l top) are consistent with the lack of mo- 
tion found by Beckers ( 1977) considering the uncertainties of 
the absolute scale of velocities (§0. 

The finding of upward motions in the inner penum- 
bra and downward motions in the out er penumbra is 
now common in the recent literature (Rimmele 1995; 
Schlichenmaier & Schmidt 2000 tldel Toro Iniesta etaLl 2001: 
Bello t Rubio et al.H2003t iTritschler et al.ll2004l) . To the best 
of our knowledge, we are the first to infer ubiquitous down- 
ward motions in a penumbra. (Recall, ubiquitous upward 
and downward motions.) This seemingly inconsistency of 
our finding with previous works is probably an artifact pro- 
duced by the spatially unresolved structure of the flows. Most 
techniques employed so far assume uniform velocities in the 
resolution element. However, observations show with insis- 
tence the presence of sever al velo city compon e nts (IB umba 
1960; Stellma cher & WiehrlfT97Tl IWjeh1fT995t iRiiedi et al.l 
1999). In our case we infer two very different types of 



motion 6 . When the observations are interpreted as a sin- 
gle resolved component, the measured velocity corresponds 
to some kind of ill-defined mean value of the actual veloc- 
ities 7 . Since (a) the true velocities can be extremely high 
(e.g., Fig.[l6|>, (b) there are both upward and downward mo- 
tions, and (c) the residual velocities are very small, measur- 
ing upflows or downflows depends on subtleties of how the 
various motions contribute to the mean values. The contri- 
bution depends on the physical properties of the penumbral 
atmospheres, as well as on method employed to measure (the 
spectral line, the portion of line profile that is employed, the 
method used to assign Doppler shifts, and so on). In order 
to show that this interpretation of the systematic flows de- 
duced by others is consistent with our inversions, we have also 

6 The finding of two types is an idealization that we impose, and it should 
be understood as the need for more than one mean component to describe the 
observed line shapes. 

7 Like the well-known convective blueshift of the photospheric lines (e.g., 
Dravins 1982). If this convective blueshift is interpreted as a net motion, one 
infers a significant global expansion solar surface, able to double the solar 
radius in less than a month. 



measured the Evershed flow of our sunspot with one of the 
techniques used in single component inferences. Downflows 
at the penumbral rim have been observed using the displace- 
ment of the Stokes I bisecto r at the line wings (e.g., Rimmele 
1995; Tritschler et al. 2004). We compute the bisector corre- 
sponding to each observed Stokes / profile of Fe I A6301 .5 A, 
the line with lowest Zeeman sensitivity in the spectra. The 
azimuthally averaged shifts of the bisectors are represented 
in Fig. [T7b (the solid line with asterisks). (The azimuthal 
average of the line-of- sight velocity cancels the contribution 
of the horizont al motions, so that only vertical v elocities re- 
main; see, e.g., Schlichenm aier & Sch midt 2000.) The simu- 
lated bisector-based measurement shows upflows in the inner 
penumbra and downflows in the outer penumbra, in agree- 
ment with single component estimates (several representative 
measurements from the literature are included in Fig. [171. The 
bisector-based mean velocity is qualitatively different from 
the true mean velocity of our inversions (the solid lines with 
error bars in Fig. H7I . This OCF weighted mean velocity 
shows downflows everywhere, meaning that the minor com- 
ponent dominates the average. However, the nature of the 
curve changes when the mean is computed weighting each 
component with its mass (i.e., mass density times OCF). Then 
the vertical velocity is not far from to some of the mean veloc- 
ities of the literature (compare the plain solid line in Fig.fTTb 
with, e.g., the dotted line). One may think that the scatter 
among individual velocities in Fig. [16| questions the small 
mean velocities shown in Fig. El However, small differences 
of mean velocities are significant since hundreds of indepen- 
dent inversions contribute to each mean value. The scatter of 
the mean decreases as the square root of the number of pix- 
els contributing to the mean, a factor of the order of 30 at 
the penumbral edge. The plain mean velocities in Fig.fTTlin- 
clude these formal error bars, supporting the reliability of our 
mean velocities. (The errors for the mass weighted mean are 
similar.) Formal error bars for the mean velocities are also 
included in Fig.[T6l 

4.6. Net Circular Polarization 

The presence and behavior of Net Circular Polarization 
(NCP) created by individual spectral lines has stimulated 
some of the important advances on the penumbral magnetic 
field structure. Any consistent penumbral model should be 
able to reproduce its properties. Our semi-empirical sunspot 
does it. The NCP is defined as 



NCP= / V(A)dA, 



(8) 
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Fig. 15. — Horizontal U r and vertical U z velocities of the two magnetic components at logr c = -1. The equivalence between the color code and the velocity 
differs for the different plots, as it is indicated by the vertical bars. Positive vertical velocities correspond to upward motions. Positive radial velocities correspond 
to outward motions. First, note the large velocities of the minor component (the two figures at the bottom), as compared to the velocities of the major component 
(the two figures at the top). Second, note the downward motions of the minor component plasma occurring over the entire penumbra. The minor component 
signals have been set to zero in the umbra, since the minor component does not contribute to the umbral spectrum. The white contour outlines the border of the 
continuum sunspot, i.e., the transition between the penumbra and the surrounding photosphere (c. 



umbral spe 
.f. Fig.CJ. 



with the integral extending a full line profile included between 
wavelengths Ai and A2. Figure[HD shows a scatter plot of syn- 
thetic versus observed NCP of Fe I A6302.5 A. Although the 
observed values tend to be larger than the synthetic ones, the 
syntheses grasp the observed behavior. The maximum value 
is around -3 mA and the signs are preserved. In addition, both 
the observed NCP and the synthetic NC P follow the ge neral 
rules found bv piling et all (Il974albl) andlMakital (119861) . and 
summarized in Sanchez Almeida & Lites ( 1992, § 4); the ex- 
treme NCP is negative (as the sunspot polarity) and occurs at 
the limb- side penumbra where the Stokes V profiles show the 
cross-over effect (profiles with three or more lobes like those 
in Fig.|3ji). 

In order to add one more piece of evidence support- 
ing the consistency of the model sunspot, we synthe- 
sized the NCP produced by Fe I A 15648 A. The NCP 
of this line is peculiar since it does not obey the rules 
of the NCP mentioned in the previous paragraph, in par- 
ticular, the NCP is not symmetric with respect to the 
line dividing the sunspot along a solar radial direction 
(iSchlichenmaier & Coll ados] EflQa Fig. 8, § 4.3). Due to 
this symmetry breaking, the extreme NCP for Fe I A 15648 A 
does not occur in the penumbral regions closest to the so- 



lar lim b. T his behavior i s attribu ted by|Schlichenmaier ~t alJ 
(2002) and iMuller et al. (2002) to the action of Magneto- 
Optical Effects (MOE) in the radiative transfer. MOE 
are significant only for very large Zeeman splittin gs (e.g., 
Landi Degl'Innocenti & Landi Degl'Innocenti 1973), which 
is indeed the regime to be applied to the strongly Zeeman sen- 
sitive Fe I A 15648 A line. Due to the existence gradients of 
magnetic field azimuth along the line- of- sight, MOE create 
NCP out of linear polari zation signals dSchlichenmaie r'etal] 
2002; Miiller et al. 2002). These gradients of azimuth are pro- 
duced by the existence of two magnetic components in the 
resolution element having different magnetic field inclination 
with respect to the vertical. When the line-of- sight is not 
along the vertical direction, the differences of vertical incli- 
nation render differences of azimuth in the plane perpendicu- 
lar to the line-of- sight. (For a very in tuitive representation of 
the effect, see Fig. 4 in | Kalmanl 19911) Th e magne tic geome- 
try re quired by Schlichenma ier et alJ (120021) and Miiller etall 
(2002) is already present in the model sunspot that we re- 
trieve and, in principle, it has the potential to reproduce the 
observed Fe I A15648 A NCP. Figure [19| shows the NCP pro- 
duced by Fe I A 15648 A assuming an axi-symmetric sunspot 
located at /i = 0.7. (We have used one randomly chosen model 
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Fig. 16. — Horizontal U r and vertical U z velocities of the two magnetic 
components of the model MISMAs at logr c = -1. We represent the mean 
(lines) and the standard deviation (shaded regions) for all pixels at a distance 
from the sunspot center. The major component is in the upper plot whereas 
the bottom plot corresponds to the minor component. Positive velocities are 
upward (U z ) and outward (U r ). The vertical bars correspond to the formal 
error bars of the mean velocities. 



MISMA to represent the thermodynamic state and the veloc- 
ities in the atmosphere of this axi-symmetric sunspot.) The 
synthetic Fe I A 15648 A NCP turns out to be in good agree- 
ment with the observed spatial distribution. Figure [19] also 
includes the NCP of Fe I A6302.5 A, which is symmetric with 
respect to the line along the solar radial direction (arrows in 
Fig.[T3. 

4.7. Wilson depression 

The presence of a strong magnetic field evacuates and 
cools down the sunspot atmosphere with respect to the non- 
magnetic quiet Sun. The resulting drop of density produces 
a depression of the observed sunspot photosphere with re- 
spect to the non-magnetic photosphere. Figure |20| shows the 
depth of the layer r c = 1 as a function of the pos ition in the 
sunspot (i.e., the Wilson depression). Refer to $ 13.51 for the 
estimate of absolute heights. The Wilson depression of the 
sunspot core is smaller than the values reported in the litera- 
ture, even keeping in mind the 100 km uncertainty worked out 
in Appendix IaI For exampl e, Allen's Astr ophysical Quanti- 
ties quotes (-0.6±0.2) Mm jSolankill2000L whereas we get 
-0.25 Mm instead (Fig. l20b). Despite the discrepancy, we 
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Fig. 17. — Average horizontal (a) and vertical (b) velocities as retrieved 
from the MISMA inversions. Positive means radially outward (a) and verti- 
cally upward (b). The solid lines with error bars show the volume average of 
the two components, whereas the plain solid lines correspond to the density 
weighted average. The other lines are shown for reference, and the y corre- 
spond to works finding down flows within the sunspot bo undary I Rimmele 
1995, the dotted dashed line; Schlichenmaier & Schmidl 2000, the dashed 
line, where we only represent one among the four targets; Tritschler et al. 
2004, the dotted line). The solid line with asterisks shows the mean vertical 
velocities we deduce from the bisector of Fe I A6301.5 A, simulating one of 
the measurements from the literature. The radial distance r has been referred 
to the sunspot radius R$ to allow comparison between sunspots of different 
radii. 



trust our estimate because the semi-empirical model sunspot 
reproduces the observations on which the large Wilson de- 
pressions reported in the literature are based on. 

First, there are estimates based on the correlation between 
magnetic field strength B and continuum intensity I r (e.g., 
IKodd & Rabinll 99llMartinez Pillet & VazQue7fl99l: 

The symbol P(Wd) stands for the gas pressure at the Wilson 
depression in the non-magnetic Sun, where I c = I c q. The re- 
lationship (|9l> is inspired by the condition of magnetostatic 
equilibrium of an axi-symmetric sunspot (for details, see the 
papers cited above). The Wilson depression Wd is determined 
carrying out a linear regression B 2 versus I c that provides a 
slope m and a constant k, 

I = mB 2 + k. (10) 
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Fig. 18.— Net Circular Polarization (NCP) of Fe I A6302.5 A. Scatter 
plot of synthetic NCP versus observed NCP. The synthetic values slightly 
underestimate the observed ones (the solid line corresponds to equal synthetic 
and observed NCPs). However, observations and modeling follow the same 
tendency and have the same spatial distribution across the sunspot (see main 
text). 
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Fig. 19.— Net Circular Polarization (NCP) of Fe I A 15648 A (left) and 
Fe I A6302.5 A (right) when one of the model MISMAs is used to produce 
an axi- symmetric magnetic field distribution. This axi- symmetric distribu- 
tion is observed from the side to mimic a sunspot observed at fi = 0.7. Note 
how the NCP of Fe I A6302.5 A is almost symmetric with respect to the 
line along the solar radial direction (indicated by the arrow, whose tip points 
the solar limb direction). This symmetry disappears for Fe I A 15648 A, 
an effect observed by Schlichenmaier & Collados 1 2002), and explained by 
Schlichenmaier et al. < 2002 ) and Muller et al. < 2002 ) in terms of a thick hori- 
zontal fluxtube embedded in a more vertical magnetic background. Our mod- 
els also explain this observation. See text for details. The spatial coordinates 
are given in arbitrary units. 
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Fig. 20. — Wilson depression as a function of the radial distance within 
the sunspot. The shaded region corresponds to the standard deviation among 
those points having the same radial distance from the sunspot center, (a) Wil- 
son depression for the observed sunspot, with fi ~ 1 . (b) Wilson depression 
if the sunspot were close to the solar limb (/i = 0.1). Note how the Wilson 
depression increases toward the limb. 
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Then, inverting equation 

W D = p- l [-k/{Aixm)], (11) 

where the depth corresponding to a give n gas pressure (the 
inverse P function or P~ l ) is derived from lSpruitl (119741) con- 
vection zone model. Figure|2]shows the scatter plot between 
the observed continuum intensity and the mean magnetic field 
strength provided by the inversion at an optical depth repre- 
sentative of the formation of the spectral lines (logr c = -1). 
It also shows two linear fits obtained using I c < 0.8 (the solid 
line) and I c < 0.5 (the dashed line). The application of equa- 
tion dTTl renders Wp — -0.5 Mm in both cases. This figure is 
consistent with the values reported in the literature, but over- 
estimates the true Wilson depression of the model sunspot by 
a factor two. 

Second, most Wilson depression measurements have been 



Fig. 21. — Scatter plot of the square of the magnetic field strength B 2 
versus continuum intensity. Linear fits of this empirical relationship have 
been used to estimate the Wilson depression of sunspots at the solar disk 
center. The two linear regressions consider continuum intensities < 0.8 (the 
solid line) and continuum intensities < 0.5 (the dashed line). Both provide an 
apparent Wilson depression of about -0.5 Mm, which overestimates the true 
Wilson depression of the model sunspot (-0.25 Mm). 



carried out for sunspots at the solar limb. Close to the limb the 
Wilson effect becomes conspicuous, and it allows estimating 
the depression from simple geometrical arguments. In partic- 
ular, the center- side penumbra becomes shorte r than the limb- 
side penumbra (e.g., Bray & Loughhead 1974, §3.8). The es- 
timates of the Wilson depression based on this foreshortening 
also exceed the value that we infer at disk center. However, 
the Wilson depression of our model sunspot varies with the 
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position on the disk. It increases toward the solar limb. The 
physical basis is simple: pressure and density scale heights 
are proportional to the temperature and, therefore, they are 
much smaller in the umbra than in the non-magnetic photo- 
sphere. Consequently, the range of atmospheric heights con- 
tributing to the observed light is narrower in the umbra. The 
umbral photosphere is not only deeper but thinner. Since limb 
observations reveal high photospheric layers, the Wilson de- 
pression necessarily increases toward the limb. We have es- 
timated the center-to-limb variation of the Wilson depression 
as follows. For plane parallel atmospheres, the optical depth 
of an atmospheric layer changes with the position on the disk 
as 

r c (/i) = r c (l)//i, (12) 

with 11 the cosine of the heliocentric angle 0, and r c (l) the 
optical depth at disk center. The Wilson depression at fi is 
the height in the atmosphere having r c (/i) = 1, which cor- 
responds to the disk center optical depth equals fi. In gen- 
eral, the plane parallel approximation is not valid to describe 
our model sunspot close to the limb, since the photon mean 
free path becomes larger than the pixel size at the disk cen- 
ter. However, it is a fair approximation to describe the Wilson 
depression of the umbral core, where the physical properties 
vary on a scale much larger than the resolution element (e.g., 
the Wilson depression is rather uniform in the central 5 Mm, 
see Fig. l20h). The same approximation is also valid in the 
quiet Sun, which we hav e us ed to set the absolute height of 
the Wilson depression (§ 13.5b . Figure l20b shows the Wilson 
depression at /i = 0.1 (0 ~ 84°) following from equation (Tl2l> . 
It turns out to be larger than the depression at the disk center, 
explicitly, close to -0.4 Mm in the umbral core. Is this value 
compatible with the observations of the Wilson effect near 
the limb? It provides a fair description of the most conspic- 
uous property of the Wilson effect, namely, the differential 
foreshortening of the limb- side penumbra d\ as referred to the 
center-side penumbra d c (see Fig. |22|>. Using simple geomet- 
rical arguments, it is possible to write down a expression for 
this foreshortening that depends only on the depression of the 
umbra Wd and the radii of the umbra Rjj and the penumbra 
Rs, 
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(Rs-Ru)^ + W D ^- 
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(Expressions similar to equation IT3l are used to in fer Wilson 
depre ssions from observed foreshortenings, e.g., Chistyakov 
Il962h Using equation fT3l and the Wilson depressions of the 
model umbra at various heliocentric angles, we compute the 
differential foreshortening. It is shown in Fig. |23] The same 
figure includes values report ed in the literature (the square 
symbols labeled as observed; Jensen et al. I ll 9691) . The model 
sunspot follows the observed tendency. The ratio is a little 
too small, but keep in mind that the observed curve represents 
mean values, and the individual sunspots do show significant 
deviations from these mean values (see., e.g., the points in 
Fig. 6 of Jensen et al. 1969). If one fits the foreshortening 
observed at the limb using as a constant free parameter Wd in 
equation JT3l . Wd — -0.45 Mm (see the dashed line in Fig. 
l23l . This value is close to that in Fig.l20b. which allows us 
to argue that the low Wilson depression of our model sunspot 
at /i = 1 is consistent with the Wilson depressions reported for 
sunspots at the limb (/i <C 1). 

Note yet another property of the Wilson depression in our 
model sunspot. The outer penumbra is depressed with re- 



d c y/ 





V 

/ 






_i_ 









v(3) = l 



2R n 
-2R, 



Fig. 22. — Schematic of a sunspot observed away from the solar disk center. 
The solid line represents the height in the atmosphere where the continuum 
optical depth equals one at the heliocentric angle of observation 6. The figure 
also shows various symbols used to describe the Wilson effect of our model 
sunspot, in particular, the different foreshortening of the limb-side penumbra 
di and the center-side penumbra d c . 
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Fig. 23. — Ratio between the size of the limb-side penumbra di and the 
center-side penumbra d c versus heliocentric angle of the sunspot. The crude 
estimate from our model sunspot (the solid line) provides a fair description of 
the mean observation (the symbols). The dashed line is a best fit description 
of the observations at the limb, and it yields a Wilson depression of -0.45 Mm. 
Heliocentric angles are given in degrees. 



spect to the quiet photosphere by about 100 km. The curves 
in Fig. |20| do not show a sharp transition at the penumbral 
boundary (distance ~ 17 Mm) because they represent the Wil- 
son depression of the magnetized fraction of the atmosphere. 
However, outside the penumbral border the stray light con- 
tamination of the model atmospheres occupies most of the 
resolution elements, and the Wilson depression of this stray 
light is basically zero. Consequently, there is a significant 
discontinuity of the observed layers at the penumbral border, 
which is bound to complicate the interpretation of all observ- 
ables. 

To sum up, the model sunspot obtained from the inversion 
provides large apparent Wilson depressions consistent with 
those reported in the literature, despite the fact that the true 
Wilson depression is rather small. In addition, there is a sig- 
nificant discontinuity of the depression at the penumbral bor- 
der. 

5. MAGNETIC FLUX CONSERVATION 

One can reproduce the asymmetries of the Stokes pro- 
files observed in sunspots assuming that the magnetic field 
varies monotonically with height ( Sanchez Almeidaj&LkeJ 



MISMA inversion of a full Sunspot 



15 



EM |Westendorp Plaza et aljIMl Nathew et al.||2Q0l. It 
requires atmospheres having vertical gradients of magnetic 
field inclination. However, the gradients are too large to 
fit in the global topology of t h e sunspot magnetic field 
(Sanchez Alm eida & Liteslil992l) . ISolanki et all (11993) for- 
mulate the problem pointing out that the divergence of the 
resulting magnetic field differs from zero and, consequently, 
the mag netic flux is not conserved. A way out was pro- 
posed by Solanki & Montavon ( 1993). The gradient of mag- 
netic field inclination does not have to be monotonic. A hor- 
izontal magnetic field embedded in a more inclined back- 
ground field also provides t he observed Stokes asymmetries. 
Solanki & Montavon ( 1993) use a single thick horizontal flux- 
tube to explain the idea, but a collection of thin fluxtubes 
does the same effect. The latter is ac tually the MISMA model 
proposed by Sanchez Al meid a et al. ( 1996), and the approach 
followed in this paper (§ 13. 111 . Since the conservation of mag- 
netic flux is central to assess the consistency of the retrieved 
model atmospheres, we have checked it with care. We use 
the differential formulation VB = 0, with B the magnetic field 
vector. 

The inversions do not provide details on the micro- structure 
that give rise to the asymmetries. They are optically- 
thin and the spectra are only sensitive to mean quantities 
(Sanchez Alm eida et alll996h . However, the condition VB = 
is also satisfied by the volume-averaged magnetic field, 
since integration and derivati on are linear operators and com- 
mute (Sanchez Alm eidal 19981) . Specifically, 

V<B>=0, (14) 

where the angle brackets < > are employed as in previous 
sections to represent volume average. These volume aver- 
ages are precisely the q uantities provided by the MISMA IC 
(Sanchez Almeida 1997). In order to carry out tests, we em- 
ploy equation (TT4l integrated azimuthally. Using cylindrical 
coordinates, equation (THl becomes 

^<X> + ~(r<l£>)=0, (15) 
oz r or y ' 

with B z and B r the vertical and radial components of the mag- 
netic field, and with the hat- symbol meaning azimuthal aver- 
age, 

- - 1 f 2 * 

<B Z > = —J^ <B z >d X - (16) 

Figurel24h shows the two terms of equation (IT5l. We compute 
the gradients at equal continuum optical depths (logr c = -1). 
Ideally, one should estimate horizontal gradients using mag- 
netic fields at a fixed geometrical height. However, we do it at 
constant optical depth keeping in mind that V < B >= must 
hold locally and that the pixel-to-pixel variation of the Wilson 
depression is negligible (see Appendix lAl. Our approach has 
the advantage of avoiding the noise introduced when setting 
a common scale of heights. The sign of the vertical gradi- 
ent has been switched in Fig.|24|so that the magnetic flux is 
conserved when the two curves are identical. The two terms 
are indeed equal within a fraction of G km -1 . For the sake of 
comparison, Fig. l24b also shows the same two terms corre- 
sponding to an inversion where the asymmetries are obtained 
by monotonic variation of the magnetic field inclination along 
the line-of- sight. We derive the horizontal and ve rtical gradi- 
ents from Fig. 9 of Westendorp Plaza et al. (2001). The verti- 
cal gradients are computed using the magnetic field strengths 
and inclinations at two extreme optical depths (log r c = and 



logr c = -2.8) and assuming dr c /dz = -140 km. The hori- 
zontal gradients are computed from the field strength and the 
inclination at logr c = -1.5. Figure l24b shows another inver- 
sion that reproduces the Stokes asymmetries of Fe I A 15648 A 
resorting to monotonic variations of magnetic field inclina- 
tion a long the line-of-sight. We use Fig. 15 of iMathew et alJ 
(2003), and the above dr c /dz. The MISMA inversion fulfills 
the condition of magnetic flux conservation in the penumbra 
better than the others. For the penumbral points, the mean 
absolute values of the residuals are, 

f 0.2 G km" 1 this work, 
|V < B > | - <^ 1.2 Gkm- 1 SIR, (17) 
1 4.3 G km" 1 SPINOR, 

w here SIR and SPINQR refer to th e na me of the IC employe d 
bvlWestendor o Plaza etafl (1200 lh and Mathe w et ail £003), 
respectively. 

The reason why MISMA inversions conform with the mag- 
netic flux conservation is simple. The gradients of magnetic 
field inclination required to reproduce the Stokes V asymme- 
tries are provided by the micro-structure, rather than from 
systematic variations of < B z >. There is no way two sep- 
arate small and large spatial scales when only smooth gradi- 
ents along the line-of-sight are permitted. Exactly for this rea- 
son, the same self-consistency found for MISMAs is expected 
in other inversions that allow for unresolved structure, even 
though it may not necessarily be optically-thin (e.g., those 
with a the single fluxtube embedded in more vertical back- 
ground; Bellot Rubi ol2004lBorrero et al .12 004) 

6. MASS CONSERVATION 

We have all the physical variables required to test the model 
sunspot for consistency with mass conservation. The issue 
of mass conservation, in particular at the penumbral bound- 
ary, has been traditionally regarded as one of the difficulties - 
troubles-mysteries of the Evershed effect (e . g . iThomasI 1 9941) . 
Figure |25] shows the mass per unit time and surface traveling 
upward (apU z ; Fig. l25b) and outward (apU r ; Fig. l25b). It 
includes the two MISMA components separately, as well as 
net mass flux carried together. First, note that the large ve- 
locities of the minor component actually transport little mass, 
i.e., similar to the mass transported by the very modest veloci- 
ties of the major component. This moderate mass flow results 
from the very low density of the minor component (Tig. [TO 
and, to a less extent, to its reduced OCF (Fig.[8|>. The mean 
mass transported upward by the major component is roughly 
balanced by the mean mass transported downward by the mi- 
nor component. Figure l25b shows the horizontal mass flow. 
It is positive for both the minor and the major components, 
meaning that it is directed outward. This horizontal flow de- 
creases to zero at the penumbral boundary (see the dotted line 
in Fig.l25b at a distance of 16 Mm). Is the mass disappear- 
ing? The horizontal mass flux seems to turn into downflows at 
the penumbral boundary. Note how the vertical mass flow be- 
comes negative at the external edge (the dotted line in Fig.l25h 
at 16 Mm). More specifically, the upflows associated with the 
major component turns to zero so that only the minor compo- 
nent downflows remain (the dashed line in Fig.l25h). Accord- 
ing to this picture, at the penumbral boundary the major com- 
ponents cease to carry mass. We just see a return mass flow 
similar to that existing in the rest of the penumbra. Formal er- 
rors for the mean net mass flows are inc lude d in F ig.l251 They 
are computed foll owing the recipe in § 14. 31 and § 14.51 As we 
point out in § 14.71 the Wilson depression changes abruptly at 
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Fig. 24. — (a) Horizontal and vertical gradients of the horizontal and ver- 
tical components of the retrieved magnetic field vector. The solid and dot- 
ted lines represent the horizontal and the vertical gradients, as indicated by 
the inset in (b). The error bars show the standard deviation of the gradi- 
ents among all points at a radial distance from the sunspot center (r/Rs, with 
Rs the sunspot radius). For the divergence of the magnetic field to be zero, 
the two curves should be equal. The top panel shows the curves resulting 
from our inversion . For comparison, the panel (b) contains the inversion by 
Westendorp Plaza et al. 1 2001), where gradie nts occur only al ong the line- 
of-sight. Panel (c) shows the inversion by Mathew et al. (2003), which have 
been carried out under the same assumption. Except for a shift, the ordinates 
of the three plots are identical. 



the penumbral border. This caveat must be considered when 
pursuing the Evershed flow outside the sunspot. 

7. THE ROLE OF THE MINOR COMPONENT IN THE 
RADIATIVE TRANSFER 

The minor component is able to modify the spectral line 
shapes, despite the fact that it carries a minute fraction of 
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FIG. 25. — Mass flows in the vertical direction (a; positive upward) and 
in the horizontal direction (b; positive outward). Despite they have very dif- 
ferent velocities, major and minor components transport similar amounts of 
mass. Note that there is no net vertical mass flow, except at the penumbral 
border (a; the dotted line). The shaded areas represent the standard deviation 
of the minor and major mass flows. The vertical bars represent the formal 
error bars of the Minor+Major mean curves (the dotted lines). 



the major component mass (about 5%; § 14. 51) . This appar- 
ent contradiction admits a reasonable explanation. Whenever 
the minor component contends with the major component to 
modify a line shape, it loses. This is the reason why the minor 
component cannot modify the cores of the lines, whose ab- 
sorption is provided by the major component. However, due 
to the Doppler shift produced by the large velocities of the 
minor component, it often absorbs in the wings of the ma- 
jor component absorption, where the contribution of the ma- 
jor component is not so overwhelmingly dominant. Thus the 
minor component can change the line wings, producing line 
asymmetries from which we infer its presence. Figurel26ltries 
to illustrate this simple explanation showing how the Stokes 
/ bisector is modified by the velocity of the minor compo- 
nent. The figure shows a Fe I A6302.5 A Stokes / synthetic 
profile corresponding to a point in the limb- side penumbra. It 
also includes the bisector, i.e., a line joining the mid points 
of the profile at each intensity level. The shift of t he bisec- 
tor is often used to measure velocities (e.g., § 14.5b . Figure 
l26b shows various synthetic bisectors obtained with the same 
model MISMA when the minor component velocity is varied 
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Fig. 26. — (a) Stokes / profile and its bisector for a typical point in the 
limb-side penumbra, (b) Bisectors corresponding to the model MISMA in 
(a) when the velocities of the minor component are modified according to 
the legend in the figure. Note how the line core is insensitive to the minor 
component speed. However, the speed of the minor component produces 
significant modifications of the bisector at the line wings. 



from zero to the value retrieved from the inversion 8 . Note 
how the bisector position at the line core is basically insen- 
sitive to the minor component velocity. However, in the line 
wings the bisector traces the minor component velocity with 
a measurable amplitude (2 km s -1 ). 

8. CONCLUSIONS 

Magnetic atmospheres having optically-thin micro- 
structure account for the line polarization observed in plage, 
network and internetwork regions (see the references in §[IJ. 
This kind of atmosphere is called Micro- Structured Magnetic 
Atmosphere or MISMA. Our paper shows how MISMAs 
also account for the Sto kes profiles of sunspots, including the 
cross-over effect ($ 13.31) . 

We have carried out a systematic fit of 10130 Stokes /, Q, 
U and V profiles of Fe I A6301.5 A and Fe I A6302.5 A. They 
correspond to the polarization emerging from a field-of-view 
containing a round sunspot (Fig.[lJ. The fits (also called in- 
versions) provide a model MISMA per pixel. Each model 
MISMA has two magnetic components with different mag- 
ne tic fi elds and mass flows interleaved along the line-of- sight 9 
(§ 13.11) . Due to the nature of the radiative transfer through 
a system of optically-thin micro- structures, such simplified 
model MISMAs only provide the mean properties of the solar 
atmosphere. In particular, we are not sensitive to the size of 
the individual micro- structures forming each component (but 
the micro- structures are optically-thin and so thinner than say 
100 km). By assembling the results of all 10130 inversions, 
we set up a model sunspot with micro- structure. Its mean 
magnetic field and mass flow is in good agreement with the 
standard sunspot magnetic field topology. The vertical mag- 
netic fields of the umbra fan out and render almost horizontal 
fields at the external penumbral border. Quasi horizontal mass 
flows of a few km s -1 exist throughput the penumbra. 

This regular large scale behavior coexists with an unex- 
pected micro- structure. We need to distinguish between the 

8 Of the order of 6 km s _1 in the layers where the line is formed. 

9 The fact that this MISMA scenario is so close to the scenarios employed 
for plage, network and inter-network regions reflects the degeneracy of the ra- 
diative transfer for MISMAs I Sanchez Almeida et al. 1996). The use of only 
two components is an idealization that we imposed, and it actually shows the 
need for more than one component to account for the observed line shapes. 



minor and the major magnetic components in each resolu- 
tion element. The major component is the one having most of 
the mass. This terminology allows us to describe the micro- 
structure of the magnetic fields and the mass flows as follows; 

1. The magnetic fields of the minor and the major com- 
ponents have very different inclinations, with the major 
component being more vertical than the minor compo- 
nent. Actually, the field lines of the minor component 
bend over to return to the sub-photosphere. This return 
of magnetic field lines occurs over t he en tire penumbra, 
rather than at the outer penumbra ($ 14.21 Fig.|9t). 

2. By construction, the magnetic field lines and mass 
flows are always parallel. Consequently, the minor 
component carries downflows. According to our inver- 
sions, both upflows and downflows are present every- 
where on the penumbra. The velocities along field lines 
are small for the major component (several hundred 
m s -1 ), and large for the minor component (10 km s -1 ; 
see Fig.[T6l». 

3. At equal optical depths, the major component has mass 
densities similar to the unmagnetized photosphere. The 
d ensiti es of the minor component are five times smaller 

(§E3I>. 

4. The minor component occupies a small fraction of the 
atmosphere (about 30 % o f the penumbra, to become 
absent in the umbra; $ 14. H and Fig. 0. 

5. Despite the large difference of velocity, the flows of the 
major and the minor components transport very similar 
masses (§|6|). In particular, there is almost no vertical 
transport of mass since the mass involved in upflows 
and downflows tends to be balanced. 

6. The minor component has magnetic field strengths 10% 
larger than the major component. If connected by field 
lines, the resulting difference of magnetic pressure can 
lead to a siphon flow from the major component to the 
minor component whos e mag nitude is similar to the ob- 
served Evershed flow ($ 14.31) . 

7. A sizeable fraction of the unsigned magnetic flux 
emerging in the penumbra returns to the photosphere 
in situ. In the external penumbral rim, the vertical mag- 
netic flux is only a small fraction of the unsigned flux 
(only 20%). This result produces a considerable differ- 
ence between th e sig ned and unsigned magnetic fluxes 
of the sunspot ($ 14.31) . 

8. The model sunspot has a Wilson depression of 
-0.25 Mm at the solar disk cent er, and it significantly 
increases toward the limb ($ 14.71) . 

Counter-intuitive as it may sound, the presence of field lines 
pointing up and down over the entire penumbra does not obvi- 
ously conflict with existing observations. We have carried out 
several tests to check the internal consistency of the retrieved 
model sunspot. 

• The MISMA inversions account for the asymmetries of 
the Stokes profiles with a magnetic field vector B which 
does not violate the VB = condition (§0). 
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• Although, the model sunspot has a Wilson depression 
smaller than the typical values found in the literature, 
it reproduces the observab les fr om which large Wilson 
depressions are inferred ($ 14.71) . 

• The model sunspot reproduces the Net Circular Polar- 
ization (NCP) of the observed lines, including their 
rules. In addition, it also reproduces the abnormal 
behavior of the NCP of Fe I A 15648 A, discovered 
and explained by Schlichenmaier and co-workers using 
magnetic structures lar ger t han the ones implicit in the 
MISMA assumption (§E5- 

• When the observed Stokes profiles are interpreted as- 
suming a single magnetic component per resolution el- 
ement, one finds upflows in the inner pen umbra and 
downflows in the outer penumbra (§ 14.51 Fig. [T7b . 
These upflows and downflows appear in the recent 
penumbral literature. According to our interpreta- 
tion, they are artificially produced w hen th e unresolved 
structure of the flows is neglected ($ 14.51) . It would be 
similar to the well known convective blueshift of the 
non-magnetic granulation, which does not correspond 
to a true systematic motion. 

• Despite of the partial cancellation of vertical magnetic 
flux in the penumbra, some net magnetic flux remains. 
The existence of such residual flux avoids an obvious 
disagreement of the model sunspot with existing ob- 
servations. Chromospheric and coronal images often 
show fieldlines connecting the penumbra with mag- 
netic structures f ar away from the sunspot (see, e.g., 
IWeiss et aljlHi- 

• The minor component contributes to the emerging spec- 
trum despite the fact that they contain a minute fraction 
of the total penumbral mass (5% or less). Such a dispro- 
portionate effect results from the large Doppler shifts 
associated with the minor components, which absorb in 



the wings of the major component absorption profiles 
(§□>• 

The interpretation of all these results in the context of a 
coherent s unspot magnetic structu re is deferred for a sepa- 
rate paper (Sanchez Almeida 2005). We prefer to disconnect 
the fact that MISMAs account for the Stokes asymmetries 
observed in sunspots (this work), from the more speculative 
physical interpretation of the retrieved micro- structure. Even 
if the physical interpretation is eventually shown to be wrong, 
MISMAs would still account for the asymmetries. Only one 
general remark concerning the global magnetic structure will 
be made. The mass flows and magnetic fluxes of the minor 
and major components are surprisingly similar. This suggests 
a physical connection between them, so that field lines and 
flows emerging as the major component return to the photo- 
sphere as the minor component of the same or a neighbor res- 
olution element. This possibility is not a direct result of the 
MISMA inversion, and it is put forward as a mere conjecture. 

The picture that emerges from our sunspot inversions is 
in qualitative agreement with the origina l MISMA sc enario 
for penumbrae put forward by Sanchez Almei daet alJ <| 1996, 
§3.1) 

A final caveat is in order. The fact that our model sunspot 
reproduces the observed Stokes asymmetries does not guaran- 
tee the accuracy of the model. However, that fact that a model 
sunspot reproduces the asymmetries should not be underval- 
ued. If the model sunspot that accounts for real sunspots is 
eventually discovered, it will be found among those in this 
rare category of models reproducing the Stokes asymmetries. 
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APPENDIX 

ERROR OF THE ABSOLUTE SCALE OF HEIGHTS 

An absolute scale of heights is derived assuming no variation of the total pressure between adjacent pixels ($ 13.51) . The purpose 
of this appendix is to estimate the error introduced by such assumption. 

The true total pressure P depends not only on the height in the atmosphere, z, but also on the horizontal coordinates. For the 
sake of simplicity, we consider an axi- symmetric sunspot so that the horizontal coordinates are reduced to the distance to the 
sunspot center r. Let us use P(r, z) to denote the distribution of pressures inferred assuming that there is no spatial variation of 
P(r,z). The estimated pressures are identical to the true pressures except for a shift of the vertical scale, 

P(r,z) = P(r,z-zbW). (Al) 

This shift, zo(r), is the error of the vertical scale that we want to estimate. The P(r,z) has been constructed so that its variation 
with the horizontal position is minimum, i.e., 

dP(r, z) = OP (r, z - zo(r)) dP (r, z - Zo(r)) d Z p(r) ^ Q 
dr dr dz dr 

This expression provides an ordinary differential equation for the error, 

dzo(r) dlnP(r,z-Zo(r))/dr 



dr dlnP(r,z-z (r))/dz : 



(A3) 



which can be integrated to evaluate zo- The right-hand- side of the previous equation can be estimated from the magnetostatic 
equilibriu m equat ion. It couples the spatial gradients of gas pressure P g , radial magnetic field B r , and vertical magnetic field B z 
re. g ..IPrie^tlfT982t §8.4), 

dr 4tt V dz dr J 

Using the conservation of magnetic flux in an axi- symmetric sunspot with radial magnetic fields, 

dB z 1 d(rB r ) 
oz r or 



and the definition of total pressure, 
it is possible to transform equation (IA4b into 



d2 d2 

P=P ' + t + t> (A6) 

dP __B 2 r _^B 2 di ^ A7) 
dr 4irr 4ty dz ' 

where i stands for the magnetic field inclination with respect to the vertical direction (Fig. [K)l>. (Note that we drop from the 
expressions the dependence of the physical variables on height and radial posi tion. We do it for the sake of simplicity, without 
compromising clarity.) The second term of the right-hand- side of equation dA7l» vanishes for the kind of atmosphere that we deal 
with (whose mean vertical gradient of inclination is zero), therefore, 

dlnP 2 sin 2 *, 

— — ~T~7> , ( A8 ) 

dr 1 + p r 

where we use (3 for the ratio of gas pressure to magnetic pressure, (3 = SirPg/B 2 . Using this equation plus the definition of the 
total pressure scale height H P , 

dz 

one can find a very compact expression for the horizontal gradient of the error, 

dzo(r) 2H P sin 2 1 

~dr~ ~ T+f3~T' (A10) 



20 



Sanchez Almeida 



100 r 




0.0 0.2 0.4 0.6 0.8 1.0 

r/R s 

Fig. All . — Error of the absolute vertical scale due to the assumption that the total pressure does not vary between neighbor pixels. They are represented 
for each radial position along the sunspot (r/Rg, with R$ the sunspot radius). The solid line is the error relative to the sunspot center, whereas the dashed line 
correspond to the error between two adjacent pixels. Heights are given in km. 



This differential equation can be easily used to estimate both, the error between adjacent pixels separated by Sr, 

S Z0 (r)^Sr^p-, (All) 
dr 

and the total error from the sunspot center to a given radial distance r, 

Azo(D= f^V (A12) 
Jo dr 

Figure lA27l shows both Szo and Azo computed according to equations (IA10L (lAllb . and JA121) . under the assumption (3=1, 
Hp = 100 km, Sr ~ 0.75", and 

* = 80°-^, (A13) 

with R s the sunspot radius. Equation dA13l» roughly describes the mean inclination of a sunspot magnetic field (e.g..lSolanki 
2000, or Fig.|9l>. The rest o f parameters are typical of our penumbral inversions. Note th at sin 2 ^/r ex r in the umbra, so that 
integrand of equation JA12I> is only significant in the penumbra. According to Fig. IA27I the pixel-to-pixel error is negligible 
(smaller than 4 km). The variation across the full spot is larger, but it does not reach 100 km. 



